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Abstract 

The joint influence of the polariton effect and Kerr-like nonlinearity on propagation 
of optical pulses is studied. The existence of different families of envelope solitary wave 
solutions in the vicinity of polariton gap is shown. The properties of solutions depend 
strongly on the carrier wave frequency. In particular, solitary waves inside and outside 
of the polariton gap exhibit different velocity and amplitude dependencies on their 
duration. 

pacs42.65.Tg, 42.25 Bs 

1 Introduction 

During the last years, due to fast progress in the fabrication of microresonators, organic and 
inorganic quantum wells etc., great interest has been attracted to investigation of electro- 
magnetic properties of these new objects, including propagation of nonlinear pulses in such 
structures (see, e.g., PP-|1]). For example, in jl] propagation of nonlinear pulses along a quan- 
tum well imbedded in a microresonator was studied for two types of nonlinearities — Kerr-like 
nonlinearity applied to envelopes of long enough pulses, and self-induced transparency (SIT) 
nonlinearity applied to short and intense pulses with frequency close the two-level resonance. 
In this study, the authors have chosen such conditions of propagation that one can neglect 
the polariton effect of formation of a gap in the dispersion law of an electromagnetic wave 
coupled with the wave of polarization in a medium. Then the problem of pulse propagation 
can be reduced to either nonlinear Schrodinger (NLS) equation, or sine-Gordon equation 
with well-known soliton solutions. However, the region of frequencies in vicinity of the po- 
lariton gap is very important, because some properties of the structures under consideration 
show up in this region only. 

One should note that related problems were already studied long time ago in the theory 
on nonlinear pulse propagation through a medium in vicinity of exciton resonances. In papers 
j3] , the polariton SIT pulses were found, but later there was made a statement [HI EI that 
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polariton effect prevents the existence of the SIT pulses. This contradiction was clarified in 
Refs. 0, where it was shown that polariton solitons do exist due to subtle balance of small 
effects, and therefore these solutions may easily be overlooked if too rude approximation 
is made in the evolution equations. Later the authors of confirmed in general this 
result and described some additional remarkable properties of the polariton SIT pulses in the 
vicinity of polariton gap beyond the perturbation theory. Similar problems have also been 
studied for the case of Kerr nonlinearity (see, e.g., Ref. and references therein). However, 
the approximations used were not justified well enough, and some properties of polariton 
solitons remained unclear. In paper [T2j, for long pulses and frequency of the carrier wave far 
enough from the polariton gap, the governing equations were reduced to the perturbed NLS 
equation for the envelope and corresponding soliton solutions are described. A closely related 
problem, namely a pulse propagation in Kerr-nonlinear medium with singular dispersion 
relation was studied in where bright and dark solitary wave solutions were found in 
vicinity of linear resonance. 

In this paper we perform a thorough investigation of polariton solitons for the case of 
Kerr nonlinearity following the method developed by Akimoto and Ikeda jH] and show the 
existence of localized solutions both inside and outside the polariton gap. 



2 Main equations 

We start with the standard equations of the classical theory of electromagnetic waves prop- 
agating through an isotropic medium (see, e.g., [Ej), 

+ uj^P + x\PfP = aE, (2) 



where eo is the background dielectric constant, and 

a = eo{ujl - UT)/4:n. (3) 

These equations describe interaction of the electromagnetic field E and the polarization wave 
P due to the Kerr-like nonlinearity (measured by the parameter x) ■ Here we ignore the effect 
of damping. The parameters, and ul characterize the dispersion law u = uj{k) of linear 
waves, which is given by the equation 

^ = 0. (4) 

As follows from (0)) the dispersion law has a gap in the frequency interval 

ojt < oj < ui, (5) 

where linear waves cannot propagate. As was mentioned above for frequency far enough 
from the polariton gap one can introduce the envelope function and the system ((H), Q 
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can be reduced to the NLS equation possessing well-known soliton solutions. Here, we are 
interested in solutions of the system (^Q), (j2)) for frequency near and inside the polariton gap 
©• 

We look for the solutions in the form of stationary linearly polarized waves, so that E 
and P may be considered as the following scalar functions 



E{x, t) = £{t - x/V)e'\ P{x, t) = [u{t - x/V) - iv{t - x/V)] e'" 
where V is velocity of the pulse, and the phase 6{x,t) is 

e{x, t) = kx-ujt- (t){t - x/V). 



(6) 



(7) 



Substitution of Eqs. (jHl), © into Eqs. (|H), ^ leads to the following system of equations for 
the variables u, f , £, 



u 



LO — ujrp -\- 2uo(j) + (j) \ u — 2[lo + (j))v — (j)v + x{u + v )u = at, 
- (u;^ -ujt + 2uj(t) + 0^) + 2{uj + (j))u + 0m + + v'^)v = 0, 



(y2 



£ 



u — uiuj + 



— (hv — 2{u + (j))v 



k epg; i 

V c2 



J_ _ a 

y2 c2 



(47r/c2 



V — viuj + 



+ + 2{uj + 0)m 



(8) 
(9) 

(10) 



where the overdot stands for the derivative with respect to = t — x/V . 



3 Linear approximation 

We suppose that at the tails of the pulse (i.e. at infinite |^|) the variables m, £, go to 
zero, so that in these regions the system (jH|)- (fTT|) can be linearized: 

ii — {uj'^ — u — 2ujv = a£, (12) 
V ~ (u^ -uj'pjv + 2uju = 0, (13) 

2 - ^) ^ = (4^/^') - + 2^ • (15) 

For the exponential dependence 

{£,u,v) = {£o,uo,vo)exp{-\^\/T), at |^| ^ oo, (16) 

where r is the duration of the pulse, the system (fT^ - (fTKjl reduces to algebraic equations 
which define the "dispersion law" and the velocity of the pulse as functions of r. It is 
convenient to introduce the variables 

X = ck/uj, Y = c/V, (17) 
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and to define the characteristic parameters 

s = I/UT, = - lu^)t^, 

so tliat the ratio 



„2 



IjJt — UJn 



does not depend on r. The variable 



(19) 



(20) 



measures the frequency in vicinity of the polariton gap. The cases of > 1 and A^ < 
correspond to upper and lower polariton branches respectively. By substituting Eq. (fTH|) into 
Eqs. (fT^ - (fTHjl . we arrive at the system 



(A2-A2)^xo- 25^2^0 = (eo/47r)£:o, 

(A2 - A2) Vq + 25^2^0 = 0, 

(X2 _ s^Y^) £, = eo(l - s^)£, + 47r[(l - s^)u, + 2sv^], 
2s {XY - eo) So = 47r[(l - s^)vo - 2suo]. 

Eliminating uq and Vq from Eqs. (j?H) gives the following system 



X2 _ ^2y2 

XY 



l_s2 _^ (l-s2)(A^-A2)-4s^n2 



(A2-A2)2+4s2f74 
(A2~A2)2+4s2Q4 



A2-A2 + (1 



the solution of which with account of Eq. (fT^ yields 

, 2 



ck 



(1-s2)(A2-A2)+4A2 
(A2-A2)2+4A202 



(A2~A2-l)2+4A2n2 
(A2-A2)2+4A2n2 



:2i) 



(22) 



(23) 



2s2 



M _ „2N , (l-s2)(A2-A2)+4A2 

^ ) ^ (A2-A2)2+4A2n2 



(A2-A2-1)2+4A202 

(A2-A2)2+4A2Q2 



(24) 



In the limit of a uniform wave (r — oo), when A^ — >■ 0, — 0, Eq. (j23|) reproduces the 
dispersion law of linear plane waves, 



ck 



eo 1 



1 

A2 



> 0, 



(25) 



which after substitution of Eq. pn|) can be transformed into the standard form. To find the 
velocity of the envelope of a linear wave, we have to take the limit ^ of Eq. (f^^ which 
gives 



V 



ep 

A6(A2 - 1) 



(26) 
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As one could expect, this velocity V coincides with the group velocity of propagating linear 
plane waves with the dispersion law (0)). 

In Eqs. (I23p and (|24p three parameters, A^, s^, and depend on the frequency u. For 
further investigation it is convenient to express and fl^ in terms of A^, 



A^ + -2 '2 



K 



and then equations ()23|) and 

{if ^ 



A2 + fi;2' 

take the form 



UJt — OJ. 



,2 ' 



1 - 



(1-A2 /{A2+k2)) (A2-A2)+4A2 
(A2-A2)2+4A2(A2+k2) 



+ + A^W 



(A2-A2-1)2+4A2(A2+k2) 
(A2-A2)2+4A2(A2+k2) 



(27) 



f28^ 



2A2 



^2 (1-A2/(A2+k2))(A2-A2)+4A2 
~ + A2+k2 (A2-A2)2+4A2(A2+k2) 



A2 



A2+k2 



(A2-A2-1)2+4A2(A2+k2) 
(A2-A2)2+4A2(A2+«;2) 



(29) 



where we have also put eo = 1 that is equivalent to the replacement c cj We remind 
that in Eqs. (|28|) and (|29|) is a constant determined by the system under consideration, 
the parameter A^ measures the frequency of the wave and the parameter A^ measures the 
duration of the pulse. 

Linear uniform waves cannot propagate with frequency within the polariton gap (0), or 



< A^ < 1. 



(30) 



However, at finite values of r, two branches of the dispersion curve are joined into one curve. 
Plots of A^ against [ckjiSf' at several values of A^ are shown in Fig. 1. As we see, these 
curves depend essentially on the values of A^, and only at A <^ 1 and far enough from 
the polariton gap we can apply the usual approach with transition to NLS equation for the 
envelope function. The velocity parameter (c/V^)^ of the pulse as a function of A^ at several 
values of A^ is shown in Fig. 2. It has real values even at frequencies inside the polariton 
gap (j30|) . It is important to note that if general nonlinear equations have a pulse solution 
of the form (jH))^© , then its velocity must coincide with its "linear approximation" ()29|1 
calculated for the tails of the pulse. 



4 Soliton solutions 

To find soliton solutions, we return to the exact equations (jS|)-(fTT|) replacing differentiation 
with respect to ^ = t — xjV by differentiation with respect to C = ^/t. Taking into account 
Eqs. (fT7j) - (p?|) . we arrive at the system 

k^u- (A2 + 2AVA2 + K^i^ + A202)m _ 2(AVA2 + /s:2 + A^^)?) 
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K^v - (A^ + 2AVA2 + + A202)^ +2(Av/A^T^ + 



where 



+A20U + + v'^)v = 0, 



A^ 



- 1 



+2Av/A2 + /t2 



- <! (A2 + /t2) 

+ A2 



cfe 



c k -| 



>£ 



A^u - m(VA2 + k2 + A0)^ 



-A2t;0 - 2(AVA2 + /t2 + A2 



0^ + 2{Av/A^T^(fg-l 



A2 
+A2 
= Air 

+K^u(P + 2(AVA2 + fi;2 + A2 



A2i} - t;( VA2 + /t2 + A0)2 



)u 



(32) 



(33) 



(34) 



(35) 



and overdot stands now for derivative with respect io C, = ^/t. 

Let us consider long pulses with A^ ^ A^. To this end, we introduce a small parameter 

e by 



A" = £"A 



2 a2 



£ < 1, 



(36) 



so that Eqs. (p?T|l - (P^ can be expanded in powers of e. Since Eqs. (jSHI) and (piljl lead to 
different series expansions in different intervals of A2, we have to consider all these cases 
separately. 



4.1 Long pulse above the polariton gap 

We begin with the case of waves with 

A2>1, 



(37) 



for which a nonlinear pulse can be presented as the envelope of propagating linear waves 
with dispersion law (|^. For long pulses with 

< 1, < A^ - 1 (38) 

the coefficients in Eqs. (jSBI), can be represented as series expansions in powers of e: 



ck 



1 



A2 



A^ - 3/t2 + 4A2/t2 
A4(A2 - 1) 



1 2 

-^ + a2£ + 



A^ 



_ A^ + 2A4/t2 + ^4 2 _ 2 
- A4(A2 - 1) ^ + • • • - 72£ + • • • , 

2«:2 4/€2(A2 + 2/€2) 3 a o ^ 



(39) 

(40) 
(41) 



We suppose that functions 0, m, v can also be represented in the form of series expansions, 
£ =6" I]^o ^"^n in the same way for 0, m, v. Analysis of equations and (jSHI)- 

(|4ip shows that for self-consistency of the procedure the series expansions of the fields must 
be as follows: 

£ =e{£o + e^£2 + ...). ^ = eiO^ + e^O^ + . . 
u = £(mo + e^M2 + • • •)? V = e{evi+ e^v^ + . . ). 

So substitution of these expansions into Eqs. yields a sequence of equations for the 

coefficients £q, £2, Oq, .... In the first approximation we obtain the relations 



Mo 



47rA2 



1 ^ 2VA2 + 

£0, Vi = Mo, 



A 



(43) 



which correspond to the plane wave solution with constant amplitude. In the next approxi- 
mation we obtain from Eqs. ()31|) and ()33|) 



-£o + ^^^^^£0 - Tefe^o = ^2 + 47r A2m2 
(72 + 1 - 4^^) £0 - (A2 + K^)a2£o 
^ (£2 + 4nA^U2 - '-^9o£o) . 



2v^ 



—%£o, 



(44) 



A2 



Combination of these two equations yields after simple transformations the equation for £0: 

X 



£0 — £0 + 

If X < 0, then it has the soliton solution 

£0 = 

where 



p3 



cosh^ 



\X\ 



a2 



and 



c 



X 

V 



A4(A2 - 1) 
- A6(A2 - 1) 



(45) 



(46) 



(47) 



(48) 



this shows that V is equal to the group velocity given by Eq. (jSEI) and does not depend on 



the duration of the pulse. Since e = A/ A = 1 / {J u"^ — r) , x = x/{^ 



£ 



\ 



327r2A6a2 1 



\X\ 



T cosh 



we obtain 



(49) 



This is the NLS type soliton solution. 

Note that the above calculations were made for arbitrary values of n. However, in ma- 
jority of applications we have {ul — ujt) ^ ^t, that is 



K 



uj\ — ujj^ 2{ujl — ujt) 



> 1. 



(50) 
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Then expansions ()39 |) - ()41|) are valid only under the condition that 

K^e^ ^ l/A{ujL - ujtYt^ < 1, (51) 

which means that the spectral width of the pulse ~ 1/r is much less than the width of the gap 
{lol —lot)- In fact, the condition (j^T|) is satisfied well enough already for 1/r ^ {ujl — ujt)/'^- 



4.2 Long pulse at = 1 

In this case a pulse cannot be described by the NLS equation for the wave packet of waves 
with wave vectors around some nonzero value. From Eqs. (j^Hj) - (j^I?j) we have the expansions 



- 1 = -1 + Vl + + . . . = -1 + ai£ + . . . , 



A2 



- 1 



(1 + ^2)3/2^ _ (2«4 + f k2 + \)e^ + 
lie - 725^ + . . . , 



'c^k 



2A ( — - 1 = 2K^e - 4(1 + 2K^)K^e^ + . . . = I3ie - p^e^ + 
Vuj / 



(52) 

(53) 
(54) 



Again we look for the solution of Eqs. (jHT|) - ()34|) in the form of series expansions which in 
this case are 

u= e^^'^{uo + £Ui + . . .), v = e^^'^vi + 

In the first approximation we have 



Mo = —£o/4:7l, Vi = 2\/\ + ilQ. 

In the next approximation Eqs. ()31|) and (jH^ give 



+ X^o = ^i/47r, (1 + n'^y^'^So + Si - aiSo = -A-nui. 



Hence, 



X 



81 (1 + K^f/'^S, - So) = -{Si + 4^0, 



and we arrive at the equation for Sq: 

So = So + 



X 



s' 



167r2(l + /t2)i/2 



Thus, we obtain the soliton solution 

S 

where V is given by 



cosh 





«-f); 







327r2(l + /t2)l/2 



Ixl 



:i + /s:2)3/2 



(55) 



(56) 

(57) 
(58) 

(59) 



(60) 



(61) 



Taking into account that e = A 
terms of the physical parameters 



1/ (S 



r ) , we can rewrite Eqs. fjUUI) , (jHT|) in 



\ Ixk cosh [i (t 



^ =1 + 



ujIt 



(62) 



Thus the velocity of the pulse depends on r (curves in Fig. 2 intersect the straight line 
= 1 at different points depending on r). Although the parameter k disappeared from 
Eq. (jn2I), the expansions ((321) -(EH) valid for k 1 provided the inequality (^e)^ <^ 1 is 
satisfied. The ratio of the amplitude of the solution (jH^ to that of the solution (jiUI) is of 
order of magnitude 



[KB] 



^1/2 



that is the amplitude at the boundary of the gap is much greater than the amplitude of the 
soliton solution far enough from the gap. This means that the pulse should be intensive 
enough to deform the gap to such extent that propagation of the wave with frequency 
= 1 becomes possible. Beyond the gap there are linear waves which can propagate with 
arbitrarily small amplitudes and nonlinear effects must only compensate dispersive spreading 
of the wave packet built from linear waves. 



4.3 Long pulse inside the polariton gap 

For frequencies inside the polariton gap, 

< A^ < 1, 

we have the series expansions 



^ -1 



-1 + 



UJ 



A4(A2-l)(A2 + /t2 



-1 + + 



2A 



1 



(1-A^)(A2+k^) 
A2 



A4(l-A2) 



2k\ 4/t2(A2 + 2/t2;^3 



A3 



A7 



£^ + ...=70 + J2£ + 

e'' + ... = Pie- p3€^ + 



where 



e< (1- A^ 



e « A^ 



(63) 

(64) 

(65) 
(66) 
(67) 



Because 70 7^ 0, the soliton solution is obtained in the first approximation, so that S and 
u do not have a small factor proportional to some power of e. The equations of the first 
approximation read 

- A^M + xu^ = (l/47r)^, (68) 
1 - A2 .. 

+ £ = -4™. (69) 
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If we were able to express u in terms of £ from Eq. (|68|) and substitute the result into Eq. 
then we would obtain the equation for S having solitary wave solutions. Unfortunately, that 
can be done only numerically except for the case when 

\x\u^ < A^. (70) 

In this limit we have ^ 

- ~4^^ ~ A2(47rA2)3^' ^^^^ 

and Eq. (jU^ takes the form 
It has the soliton solution 

327r2A6(l- A2)(cji-cuf,) 1 



(73) 



\ \x\ cosh [7 (t - f 

where the pulse velocity is given by 

(^"^y-l I (l-A^)(A^ + ^^)(^^-^lK^ _ (74) 

The condition (f7n|) can be transformed with the use of the estimate u ~ a/(47rA^) into 

1 - A^ < 1, (75) 

that is the frequency must be close enough to the upper limit of the polariton gap. 
The solution (|7Hjl only applies to long enough pulses 

< 1 - Al (76) 

Note that in this case the amplitude does not depend on the pulse width r, but there is a 
strong dependence of its velocity on r. If k 3> 1, then the ratio of the amplitude (f75|) to the 
amplitude (jiUI) is of order of magnitude 

~ {ney^ ~ {lul — u)t)'t > 1, 

i.e., it is much greater, as one could expect, than the amplitude of the soliton solution 
at the boundary of the gap. 

When the condition (f7H|) is not fulfilled, we have to solve Eqs. (jF)H|) . numerically. To 
this end, let us introduce new variables 



£ = El^\xl u = Ul^\xl (77) 
where x = "Ixl^ so that Eqs. ()68l69p reduce to 

I - /S? ■■ 

A2[/ + [/3 = -(l/47r)E, —-—E + E = -A'kU. (78) 
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Numerical solution of the first equation (fTHj) yields the function 



U = U{E) (79) 

the substitution of which into the second equation of (fTHj) gives the differential equation for 
E: 

1 _ A2 ^2 p 

--^-j^ = -E-4nU{E) (80) 



A2 dC 



or 



d^E _ dV{E) 



where 




V{E) = (^^E' + An U{E')dE^ 



f82l 



may be considered as a potential in which a particle moves according to the Newton equation 
(j8ip . A typical plot of the potential (j82|) is shown in Fig. 3. The maximum value of E 
(soliton's amplitude a) is determined by the point E = a where V{E) vanishes. Then the 
solution of Eq. (IHUj) with initial conditions E{0) = a, E'{0) = provides the soliton solution. 
We have done these calculations for several values of A. The results are shown in Fig. 4 



where the dependence of the pulse profile E on the scaled variable y A^/ (1 — A^) ( is plotted. 
The dependence of the amplitude on A is shown in Fig. 5. 



5 Conclusion 

In conclusion, the nonlinear optical pulse propagation in the frequency region in the vicin- 
ity of polariton gap is studied. The problem is described by the coupled set of the Maxwell 
equations for electromagnetic field and material equation for macroscopic polarization allow- 
ing for Kerr-like nonlinearity. To solve this nonlinear system of equations analytically, the 
approach (8j based on series expansion in powers of a small parameter related to the width 
of the polariton gap and pulse duration is used. Different bright solitary wave solutions, 
depending on the position of carrier wave frequency with respect to the polariton gap, are 
found and their parameters are expressed in terms of material system parameters. Outside 
the polariton gap the soliton solution corresponds to well-known soliton of the NLS equation 
for the envelope of the wave packet made of plane waves. However, inside the polariton gap 
there are no plane wave solutions and the notion of their envelope looses its physical sense. 
Nonetheless, there the solitary wave solutions are possible with high enough values of elec- 
tromagnetic field strength so that local value of the polariton gap diminishes in the center 
of the polariton gap due to Kerr nonlinearity. The difference in physical situations outside 
and inside the polariton gap is reflected in different dependencies of soliton's amplitude on 
the pulse duration r — the amplitude does not depend on r inside the gap, it is proportional 
to r^^/^ at the gap boundary and proportional to r^^ far enough from the gap. 

S.D. thanks Marseilles University for kind hospitality and NATO Linkage Grant PST. 
CLG. 978177 for partial support. 



11 



References 

[1] E. Burstein and C Weisbuch, (eds) Confined Electrons and Photons. New Physics and 
Applications, New York, Plenum, 1995. 

[2] A. D'Andrea et al (eds) Optics of Excitons in Confined Systems, Bristol, lOP, 1991 

[3] F. Kajzar, V.M. Agranovich and CY.-C Lec, (eds) Photoactive Organic Materials. 
Science and Applications. Kluwer, Dordrecht, 1996. 

[4] V.M. Agranovich, A.M. Kamchatnov, H. Benisty, and C. Weisbuch, Nonlinear pulse 
propagation along quantum well in a semiconductor microcavity, in Beam Shaping and 
Control with Nonlinear Optics, F. Kajzar and P. Reinisch, (eds) NewYork, Plenum, 
1998. 

[5] A. Schenzle and H. Haken, Opt. Commun. 6, 96 (1972); H. Haken and A. Schenzle, 
Z. Phys. 258, 231, (1973). 

[6] E. Hanamura, J. Phys. Soc. Jpn. 37, 1553, (1974). 

[7] M. Inoue, J. Phys. Soc. Jpn. 37, 1560, (1974). 

[8] O. Akimoto and K. Ikcda, J. Phys. A: Math. Cen. 10, 425 (1977); K. Ikeda and O. Aki- 
moto, J. Phys. A: Math. Cen. 12, 1105 (1979). 

[9] S.V. Branis, O. Martin and J.L. Birman, Phys. Rev. A43, 1549 (1991). 

[10] J. GoU and H. Haken, Phys. Rev. A18, 2241 (1978). 

[11] S.A. Moskalenko, P.I. Khadzhi and A.Kh. Rotaru, Solitons and nutation in the exciton 
region of spectrum, Kishinev, Shtiintsa, 1980; S.A. Moskalenko, Introduction into the 
theory of excitons of large density, Kishinev, Shtiintsa, 1983, (in Russian). 

[12] LB. Talanina, M.A. Collins and V.M. Agranovich, SoUd State Comm. 88, 541 (1993). 

[13] U. Peschel, T. Peschel and F. Lederer, J. Opt. Soc. Am. 14, 2994 (1997). 

[14] V.M. Agranovich and V.L. Ginzburg, Spatial Dispersion in Crystal Optics and the The- 
ory of Excitons, Berlin, Springer, 1984. 



12 



Figure captions 

Fig. 1. Dispersion relation of the carrier wave for different values of the pulse duration 
r measured by the parameter A (see Eq. (HH)). 

Fig. 2. Pulse velocity as a function of the carrier wave frequency for different values of 
the pulse duration r measured by the parameter A (see Eq. (fT^ ). 

Fig. 3. Potential as a function of the electric field amplitude for A = 0.2. 

Fig. 4. Profiles of solitary wave solutions for three different values of the carrier wave 
frequency (see Eq. fOUj) ) inside the polariton gap. 

Fig. 5. Dependence of the amplitude of the solitary wave on the carrier wave frequency 
(see Eq. (HOI)) inside the polariton gap. 
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